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Abstract 

The addition of transition metals (TM) to IIFV semiconductors radically changes their electronic, 
magnetic and structural properties. In contrast to the conventional semiconductor alloys, the lattice 
parameter in magnetic semiconductor alloys, including the ones with diluted concentration (the 
diluted magnetic semiconductors - DMS), cannot be determined uniquely from the composition. By 
using first- principles calculations, we find a direct correlation between the magnetic moment and 
the anion-TM bond lengths. We derive a simple formula that determines the lattice parameter of 
a particular magnetic semiconductor by considering both the composition and magnetic moment. 
The formula makes accurate predictions of the lattice parameter behavior of AlMnN, AlCrN, 
GaMnN, GaCrN, GaCrAs and GaMnAs alloys. This new dependence can explain some of the 
hitherto puzzling experimentally observed anomalies, as well as, stimulate other kind of theoretical 
and experimental investigations. 

PACS numbers: 75.50.Pp, 61.66.-f, 61.66.Dk, 71.20.Nr, 71.20.Be 
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In the last years, the field of spin electronics or spintronics, which aims the use of the 
spin carriers in addition to the charge for electronic devices, has emerged as the new fron- 
tier for the future integrated circuit devices technology. The suggested applications include, 
e.g., spin transistors, spin valves, spin-polarized light emitting diodes, as well as, potential 
applications in the field of quantum cryptography and quantum computing [H El [3]. To 
implement this new technology, efficient techniques for spin injection, manipulation and 
detection are necessary, making it crucial the development of materials with the needed 
properties PQ [2]. Materials under intense investigation are the diluted magnetic semicon- 
ductors (DMS), in which a conventional semiconductor is doped with TM atoms The 
great advantage of DMS is to allow an interplay between magnetic and electronic properties, 
being structurally compatible with most semiconductors epitaxially grown. Primarily, great 
attention was paid to the arsenide-based DMS. Nowadays, a large class of DMS has already 
been proposed, being extensively studied, both experimental and theoretically [5] . 

Among all properties of DMS, the structural are of fundamental importance. The greater 
part of the III-V semiconductors crystallizes in the cubic zincblende (zb) structure, only the 
nitrides, that under ambient conditions crystallize in the hexagonal wurtzite (wz) structure. 
However, successful growth of zb-AlN, GaN, and InN has been reported [2J. On the other 
hand, the structure of nitride based TM is the tetragonally distorted NaCl structure, NiAs- 
like for MnAs, and MnP-like for CrAs. For III-V films containing no more than few percent 
of TM, one expects the structure to be that of the host. Moreover, there has been an effort 
to grow these binary compounds in the zb structure using molecular beam epitaxy. Recently, 
it was reported the growth of zb-MnAs [6j, and zb-CrAs [7j. The formation of these binary 
compounds with the desirable zb structure is difficult because of inherent growth problems 
[6l E]. Thus, their lattice constant in zb structure are not well established yet, and further 
experimental investigations are still needed. 

In a pioneer work, Masek et ai, by using a density functional theory (DFT) calculation, 
related the lattice parameter behavior with the quality of the samples, by taking into account 
the incorporation of TM in the interstitial positions [9]. Afterwards, the lattice parameter 
behavior (increasing with the TM concentration) was taken as a fingerprint of the presence 
of interstitial TM, as well as, arsenic antisites in the GaMnAs samples. Furthermore, once 
the lattice constant can be precisely measured by X-ray diffraction, if one knows its relation 
with the composition of the alloy, the latter is directly determined, and in DMS the accurate 
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determination of the alloy composition is of crucial importance because the TM content 
strongly affects the Curie temperature [TDj. However, it was pointed out that clearly exists 
a phenomenological correlation between a and x, but it is quite complex and not really 
understood [llj. Moreover, it was experimentally observed that the GaMnAs is an example 
of a system which does not obey Vegard's law [12j. 

Deviations from Vegard's law are derived from many factors, as electronic and/or mag- 
netic effects, as already observed for metallic alloys [T^ [T5] . The exchange integrals in 
the Heisenberg model, describing the super exchange interactions between magnetic ions in 
DMS are strongly dependent on the lattice parameter [12]. An important property for a 
material to be useful for spintronics is its possible ferromagnetism. Therefore, it should be 
of great interest, for a better understanding of the magnetic properties of DMS, a rigorous 
study of the effect of TM impurities in the structural properties of the magnetic semiconduc- 
tors. Moreover, it should be very useful to have a quantitative relation between the lattice 
parameter and the magnetic properties, which, up to our knowledge, was not addressed until 
now for magnetic semiconductor alloys, nor for its dilute concentration. 

In this work, we present a theoretical study of structural, electronic, and magnetic prop- 
erties of zb AlMnN, AlCrN, GaMnN, GaCrN, GaMnAs and GaCrAs alloys, considering the 
whole range of composition. In particular, we focus our attention on the lattice parame- 
ter, analyze its relation to the other properties, and derive a relation between the lattice 
parameter and the average magnetic moment. Moreover, we also analyze how the behavior 
of lattice parameter is related to the electronic transition from the half-metal to the metal 
phase. 

We used the frozen-core projector-augmented wave method as implemented in the "Vi- 
enna ab-initio simulation package" (VASP-PAW code) within a spin-polarized DFT [17] . We 
adopted the generalized gradient approximation for the exchange-correlation potential in the 
version proposed by Wang and Perdew [IB]. It has been demonstrated that this method is 
able to successfully describe the physical properties of several DMSs [5l |19]. In order to 
calculate the ferromagnetic (FM) state, 16-atom supercells have been used, and varied the 
number and the position of TM atoms (substitutional) in the cation sublattice, which results 
in 16 configurations not related by symmetry. The lattice constant was optimized by total 
energy minimization for each configuration and the atomic positions were relaxed. In order 
to discuss the properties of the alloys we made averages over the configurations with the 
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same composition x. 

Fig. [T] depicts the results for the lattice parameter of AlMnN, AlCrN, GaMnN, GaCrN, 
GaMnAs and GaCrAs for the whole range of TM composition. By using the values of the 
optimized lattice of FM MnN, CrN, and MnAs we note that Vegard's law predicts values 
in complete disagreement with those observed. We observe an almost linear first increase 
in comparison to the Vegard's law followed by a fast decrease of the lattice constant with 
the Cr and Mn content. The case of GaCrAs is an exception that will be explained later. 
This behavior is in agreement with other few available theoretical results [HI 1201 EI], as well 
as, with experimental data obtained by X-ray diffraction for AlMnN, AlCrN [221 123] and 
GaMnN [211 ES] in the wz structure. Here, it is worth to point out that the results for wz 
and zb should be very similar. We note that the increase of the lattice parameter is here 
due to the substitutional TM, not the interstitial as in Ref. [9]. In the case of GaMnAs, 
it is experimentally observed an increase of the lattice constant with the Mn concentration 
for low Mn content [llj, which is in disagreement with our prediction. In this case, several 
other factors can be considered as the reason for this increasing, the interstitial atoms. As 
antisites, the presence of defects due to the growth at low temperatures, and deviations from 
stoichiometry [9llIIl[I2]. 

As Vegard's law is a consequence of the size difference between the host atoms and the 
dopant, we infer that the observed deviation must result from an electronic and/or magnetic 
effect. This kind of deviation has been studied for solid solutions, but only for metal alloys 
[T3| [TH [T3] . In order to investigate this, we analyze the nearest-neighbor distance behavior 
with the TM composition, which is depicted in Fig. [2] for the AlMnN. Although the AIN 
lattice parameter suffers an expansion when doped with a low concentration of TM, the 
bond lengths TM-N remain nearly constant below a composition of ~30% TM. Another 
characteristic that we can note is the correlation between the bond length and the total 
magnetic moment, i.e., there exists an interdependence between the structural and magnetic 
properties in those alloys. For low TM concentrations the magnetic moment remains an 
integer, an evidence that the material is half-metallic [8], which is consistent with the null 
density of states (DOS) for the minority spin-component, as shown for AlMnN in Fig. [s] 
(in order to analyze the DOS for low TM-content in more detail, we also used 64-atom 
supercells). When the TM composition increases, a transition from the half-metallic to 
metallic behavior happens, as it can be observed by the non-integer total magnetic moment. 
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and by the DOS shown in Fig. |3] 

By taking into account this correlation between bond-length and the magnetic moment, 
we estimate the contribution from magnetization to the variation of the lattice parameter 
as follows. The Vegard's law is given by: 

a(iii,TM)v(a;) = (1 - a;)aiii_v + xctm-v, (1) 

where am.v and axM-v are the lattice parameters of the binaries, with III and V elements 
of column 3A and 5A from periodic table, respectively. For a zincblende structure, we can 
rewrite the equation above by expressing ctm-v in terms of the nearest-neighbor distance 

C^TM-V, 

4 

a(III,TM)v(a;) = (1 — x)aiii_Y + X-j^dTM-Y- (2) 

For semiconductor alloys, the anion-cation bond-lengths do not converge to a single value, 
but instead remain close, throughout the composition range, to their respective values in 
the pure binary compounds [26]. This is not the case for magnetic semiconductors. Here, 
not only ciTM-v dramatically changes with the TM content, but also has a direct and ap- 
proximately linear relation with the magnetic moment as 

'^TM-v(a;) = a + (3^{x), (3) 

while the bond-length of the III-V remains approximately constant {dAi-N — 1-90 A, 
dca-N — 1-97 A, and dca-As — 2.49 A). So that by generalizing Eq. [2] 

4 

a(iii,TM)v(a;) = (1 - a;)aiii_v + ^;^(" + Pl^i^))^ (4) 
from which follows immediately 

a{iii,TM)v(a;) = (1 - a;)aiii_v + xotm-v + - /ixM-v), (5) 

where the first part is exactly the Vegard's rule, and the last term represents the magnetic 
influence on lattice parameter, with C = ^ a constant related to a magnetoelasticity 
between the TM and the anion, with unit of A//xb; and /itm-v is the magnetic moment 
per TM atom for the endpoint binary TM-V. The parameters C thus obtained were all 
about 0.1 A/fiB- Explicitly they were 0.0914, 0.1199, 0.0992, 0.1307 and 0.1401 A/fiB for 
AlMnN, AlCrN, GaMnN, GaCrN, and GaMnAs, respectively. A similar relation between 
the magnetic moment and the lattice parameter was already obtained a long time ago by 
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Shiga for the 3d transition metal alloys Fe-Co and Fe-Ni [13], and Niculescu et al. for FeCoSi 
alloys [H]. The results for the equilibrium lattice parameter and the one obtained through 
the expression above are presented in Fig. [TJ We can observe that the suggested relation 
for the lattice parameter fits very well the calculated equilibrium lattice parameters. In the 
following, we give an interpretation for the lattice parameter behavior, and consequently the 
suggested equation. 

By analyzing the electronic structure for the zb alloys, the TM d ion levels are split by 
the tetrahedral crystal field into ^2 and e. Exchange interactions further split these levels 
into spin-up and spin-down levels. The ^2 levels on the Mn atom may hybridize with the 
levels with the same symmetry the anion j»-states, and can take part in the bond formation, 
while the e levels have no other states available for significant coupling, not taking part in 
bond formation. The introduction of TM impurities on the host creates a spin-polarized 
impurity band, and the Fermi level lies on it, until the TM concentration achieves a value 
for which the alloy becomes a metal. 

Let us first analyze the case of low TM content, as shown in Fig. |3]for AlMnN. In this 
case the spin-polarized band is very narrow with almost no dispersion, meaning that the 
overlap between adjacent Mn-induced levels is small, and the hybridization is not significant, 
since in this case the level is in the middle of the gap. Thus the electrons on these levels are 
localized, which reduces the bond strength, resulting in an increase of the bond-length of the 
TM-anion. This is why the TM-anion bond-length is enlarged in comparison to the bond 
lengths of the metals MnN, MnAs and CrN, where all the outer electrons do participate 
in the bond formation, resulting in a stronger bond and reduced bond-length. Therefore, 
for low TM content the alloy maintains its semiconductor character and the bond-lengths 
remain close to the enlarged TM-anion bond-lengths. At the same time, as the bond-length 
remains constant, so does the magnetic moment, and the alloy is half-metal. Now, as the 
TM content increases there is a transition from half-metal to metal phase, and there is 
not an impurity level anymore, nor a semiconductor character for the bond-lengths, but a 
band with d character, where the TM orbitals are overlapping. This overlapping decreases 
the kinetic energy, since the d-electrons become more itinerant and less localized. With 
the delocalization and decreasing of the kinetic energy, the bonding energy increases and 
the lattice parameter decreases. At the same time occurs the transition from localized to 
delocalized state, the minority spin states becomes populated, leading to a reduction of the 
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magnetic moment and of the lattice parameter. 

Finally, in order to give support to our conclusions, it should be useful to obtain the 
lattice parameter relation for an alloy which is half metal in its whole range of composition. 
This is the case of the GaCrAs alloy. In this case according to our proposition the lattice 
parameter should follow the Vegard's rule, since the magnetic moment of GaCrAs is the 
same in the whole range of composition, and equal to 3.0 fiB- Fig. [l](f) depicts the result 
for the lattice constant behavior of this alloy. We can observe that there is no transition of 
behavior as for the other alloys, and the lattice parameter nearly follows the Vegard's law 
with a small bowing parameter of 0.04 A. Therefore corroborating our predictions. 

In summary, in this work we have analyzed the lattice parameter variation of cubic mag- 
netic semiconductors of family III-V by using first-principles calculations. We conclude that 
incorporation of TM impurities can significantly affect the lattice parameters. A large de- 
viation from Vegard's law was observed. We noticed a correlation between the magnetic 
moment and the interatomic distances between the TM and its first-neighbor, which lead 
us to conclude that there is an influence from the electronic/magnetic properties. Hopefully 
these studies on the correlation of the electronic and magnetic effects with the structural 
ones, will play an important role in the understanding of the remaining questions and stim- 
ulate other kind of theoretical and experimental investigations, allowing a more thorough 
understanding of these enigmatic materials. 

This work was supported by the Brazilian funding agencies CNPq and FAPESP (procs. 
05/52231-0 and 06/05858-0). 
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FIG. 1: Equilibrium lattice parameters (squares) as function of the TM composition for the FM 
AlMnN, AlCrN, GaMnN, GaCrN, GaMnAs and GaCrAs. The sohd line represents the fit according 
to Eq. [5| and the dashed line the lattice parameter behavior according to the Vegard's law. 
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FIG. 2: (color online) Total magnetic moment (red triangles) and average bond-length (blue 
squares) for the ferromagnetic state of AlMnN as a function of Mn composition. The dashed 
line represents the integer value 4.0/2^- 
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FIG. 3: (color online) Total spin resolved density of states of AlMnN for different Mn compositions 
(solid lines). The partial DOS of 3d electrons, e (red shaded) and t2 (blue shaded) states, of the 
Mn atoms is also shown. The vertical dashed lines denote the positions of the Fermi level. 
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